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Ground-state and finite-temperature properties of a special class of exactly solvable Ising- 
Heisenberg planar models are examined using the generalized decoration-iteration and star- 
triangle mapping transformations. The investigated spin systems exhibit an interesting quan- 
tum behaviour manifested in a remarkable geometric spin frustration, which appears notwith- 
standing the purely ferromagnetic interactions of the considered model systems. This kind of 
spin frustration originates from an easy-plane anisotropy in the XXZ Heisenberg interaction 
between nearest-neighbouring spins that favours ferromagnetic ordering of their transverse 
components, whereas their longitudinal components are aligned antiferomagnetically. 
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1 Introduction 

Quantum behaviour of low-dimensional solids belongs to the most fascinating topics to be ever 
discovered in the field of condensed matter physics. Generally, it is expected that the clearest 
manifestations of quantum phenomena should emerge in the low-dimensional antiferromagnets 
which consist of spin carriers having low quantum spin number. With regard to this, much effort 
has recently been devoted to the two-dimensional antiferromagnetic spin systems that have an 
obvious relevance in elucidating the high-Tc superconductivity of cuprates [1], quantum phase 
transitions [2], geometric spin frustration [3], etc. On the other hand, the two-dimensional quan- 
tum ferromagnets remain rather unexplored from this point of view as their macroscopic quantum 
manifestations would be, on the contrary, rather surprising. 

Nevertheless, the two-dimensional ferromagnetic spin systems with a spatial anisotropy rep- 
resent very interesting systems as they may exhibit spontaneous long-range order possibly af- 
fected by quantum effects. A strong evidence supporting this statement has been provided by 
a recent exploration of outstanding quantum antiferromagnetic phase, which was discovered in 
the mixed-bond Ising-Heisenberg models (to be further abbreviated as IHMs) despite the pure 
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ferromagnetic interactions of the considered model systems [4-6]. This remarkable finding sug- 
gests that an interesting quantum behaviour can also be found in the anisotropic ferromagnetic 
systems owing to a competition between the easy-plane XXZ Heisenberg interactions and the 
familiar Ising interactions of the easy-axis type. In the present article we shall show that this 
unusual behaviour represents a generic feature of wide class of the mixed-bond IHMs that arises 
independently of the lattice topology and hence, on close-packed lattices it leads to disordered 
ground state due to a presence of peculiar geometric spin frustration. 

This paper is organized as follows. In the next section we shall briefly describe model sys- 
tems and the method that enables to obtain their exact solution. The most interesting results are 
presented and discussed in Section|3] Finally, some concluding remarks are given in Section|4] 

2 Model systems and their exact solution 

Let us begin by considering periodic planar lattices composed of two kinds of spin sites; so- 
called Ising spins /i which interact with other spins through extremely anisotropic coupling Ji 
containing their one spatial component only and the Heisenberg spins S which interact among 
themselves via spatially anisotropic coupling J(A) containing all three spatial components. The 
parameter A allows to control a spatial anisotropy in the XXZ Heisenberg interaction, which can 
be either of easy-axis (A < 1) or easy-plane (A > 1) type. The total Hamiltonian of mixed-bond 
IHMs then includes both Ising as well as Heisenberg bonds and thus, it acquires this simple form 

n = ^jY,[A{srsj + sfs^) + s^s^] - Ji ^ ^,^Af , (1) 

(ij) (fe.O 

where S" (a = x, y, z) and are spatial components of spin-1/2 operators, the first summation 
is carried out over the nearest-neighbouring Heisenberg spin pairs and the second summation 
runs over the nearest-neighbouring pairs of the Ising and Heisenberg spins, respectively. 

In what follows we shall restrict ourselves to a special sub-class of the model Hamiltonian 
Q, namely, we will further refer just to spin systems that consists of Heisenberg spin pairs 
each interacting with either two or three nearest-neighbouring Ising spins. Several examples of 
such lattices are sketched in Fig. [J where full (empty) circles denote lattice positions of the 
Ising (Heisenberg) spins. A common feature of the lattices a) and b) is that the Heisenberg spin 
pairs are placed in between two outer Ising spins (model A), while on lattices c) and d) they 
are surrounded by three outer Ising spins (model B). In view of further manipulations, it is very 
advisable to rewrite the total Hamiltonian as a sum of bond Hamiltonians, i.e. Ti = l-Lk, 
each containing all the interaction terms of one Heisenberg pair (spin clusters on Ihs of Fig.|2} 

Hk = - J[A(S'fciS'^2 + Ski^k2) + SkiSk2} " Skhi - S^h2. (2) 

Above, hi ^ h2 = Ji(Afei + Afe2) for model A, while hi = Mpli + Alg), /i2 = ■h{fj-l2 + AL) 
for model B . On account of the commutability between different bond Hamiltonians, the partition 
function of IHMs can be partially factorized into a product of bond partition functions 



N N 

ZiHM Trfc cxp(-/37ife) = Y.X{^k, 

{m} k=l {^} k=l 



(3) 
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Fig. 2. Schematic representation of an extended: a) decoration-iteration and b) star-triangle transformation. 



where /3 — l/{kBT), fee is Bohzmann's constant, T denotes absolute temperature, the symbol 
Ti'k stands for a trace over spin degrees of freedom of fcth Heisenberg spin pair and the summation 
available configurations of the Ising spins. 
A crucial step of our approach represents calculation of the bond partition function ac- 
complished by adopting the same steps as used in Ref. [4] and its subsequent replacement by 
suitable mapping transformation. For this purpose, we shall utilize generalized versions of the 
decoration-iteration [7] and star-triangle [8] transformations originally introduced by solving pla- 
nar Ising models. It is noteworthy, however, that the mapping transformations have been later on 
remarkably generalized by Fisher [9, 10] who firstly pointed out that decorating Ising spins may 
be, in principle, substituted by arbitrary statistical systems. In our case, both the transformations 
can be schematically illustrated as in Fig.|2]and mathematically they can be expressed as follows 



Zk = 2 cxp(/3J/4) cosh(/3/ii) + 2 exp(-/3 J/4) cosh(/3JA/2) = Aexp(/3i?^^i^^2), (4) 
Zk = 2 cxp(/3J/4) cosh[/3(/ii + /i2)/2] + 2 exp(~/3 J/4) cosh[f3y/ihi - + (JA)2/2] 
= B cxp[/3Ri{fil^fil^ + Hl^nl^) + f3R2nlifil2]> (5) 
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where A and R (B, Ri, and R2) are unknown mapping parameters emerging in the decoration- 
iteration (star-triangle) transformation that are unambiguously given by a self-consistency con- 
dition of the mapping relations. Unknown mapping parameters can be calculated following the 
standard procedure [9, 10], for instance, the mapping holds for model A if and only if 

A = 2exp(/3J/4)(Fit/2)'/', f3R = 2lniVi/V2), (6) 

and the functions Vi and V2 are defined as follows 

Vi = cosh(/3Ji) + exp(-/5J/2) cosh(/3JA/2), (7) 
F2 = l + exp(-/3J/2)cosh(/3JA/2). (8) 

One the other hand, the model B can be treated within the mapping transformation (|5} when 

B = 2cMPJ/'^){ViV2Viy^\ f3Ri=^HVi/V2), /3i?2 - Wl/s/lf )■ (9) 
In addition to the functions Vi and V2 defined above we have introduced here another function 



V3 = cosh(/3Ji/2) + exp(-/3J/2) cosh(^f3^ Jf + (JA)2/2) . (10) 

At this stage, a direct substitution of the transformations and (|5} into the formula (|3} 
gained for the partition function of IHMs leads to a mapping relationship with partition functions 
of the spin-1/2 Ising models on undecorated lattices (scheme from Fig.|2j, in fact, 

Model A: Zjhm = Zim{P, R), (11) 

Model B: Zihm = B^Zim{P, i?2). (12) 

Notice that the mapping relations ™d il2\ allow simple calculation of other quantities as 
well. For easy reference, we list the most important ones: sub-lattice magnetization of the Ising 
spins to| = (A^i), sub-lattice magnetization of the Heisenberg spins mf^ = (S^-^), pair corre- 
lation functions between the nearest-neighbouring Heisenberg spins C^^ = (8^18^2)' ^uh = 
and between the nearest-neighbouring Ising and Heisenberg spins Cf^ = (Afei'S'fei)- 



3 Results and Discussion 



In this section, the most interesting results for both investigated model systems will be presented 
and discussed in detail. Although the results obtained in the preceding section hold irrespectively 
of a character of the exchange parameters J and Ji, our further analysis will be for simplicity 
restricted to the ferromagnetic systems ( J > 0, Ji > 0) only. In order to avoid confusion, the 
results obtained for both investigated model systems will be discussed separately. 



3.1 Model A 

Initially, let us look more closely on the ground-state behaviour The first-order transition line 
Ac = 1 + 2Ji/J separates regions of stability of two different ground-state phases for both 
IHMs with a symmetric decoration shown in Fig.^)-b). As could be expected, there occurs a 
s,\mp\t ferromagnetic phase (FP) characterized by saturated sub-lattice magnetization = 1 /2, 
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Fig. 3. Critical temperatures of IHMs with a sym- 
metric decoration of the square (broken lines) and 
triangular (solid lines) lattices shown in Fig.©)-b). 



Fig. 4. Critical temperatures of IHMs with an asym- 
metric decoration of the triangular (Fig. \\^) and 
orthogonal-dimer (Fig.^i) lattices. 



nih = 1/2, as well as pair correlation functions C^^ = 1/4, C^^ = 1/4 when the exchange 
anisotropy is below its critical value Ac. Within FP one also finds C^^ = clearly indicating 
classical character of this long-range ordered phase. However, there also appears a remarkable 
disordered phase (DP) without spontaneous ordering (rrii = 0, m/j — 0) when the exchange 
anisotropy exceeds its critical value Ac. It should be emphasized, nevertheless, that the Heisen- 
berg spin pairs remain strongly correlated even within DP as suggested by (7^^ = 1/4 and 
^hh ~ ~ 1/4- By contrast, the nearest-neighbouring Ising and Heisenberg spins are completely 
uncorrelated, indeed, one finds Cf^ = within DP. Thus, it can be readily understood that 
DP emerges due to a peculiar geometric spin frustration originating from an antiferromagnetic 
alignment of the longitudinal components of the Heisenberg spins that accompanies preferred 
ferromagnetic ordering of their transverse components. In agreement with this statement, it can 
be easily checked that the function |5'dp) — Yii It)* life :79(ITi) + liT))fc describes the ground 
state of DP (the first product is carried out over all Ising spins and the second one runs over all 
Heisenberg spin pairs). According to this, one may conclude that all Ising spins are free to flip in 
DP on account of the spin frustration caused by Heisenberg spin pairs, which exhibit an interest- 
ing coexistence of ferromagnetic (transverse) and antiferromagnetic (longitudinal) correlations. 

Now, let us step forward to the finite-temperature behaviour. The variation of critical temper- 
ature of IHMs with a symmetric decoration of the square (Fig.^) and triangular (Fig.fQ)) lattices 
is displayed in Fig.|3]for two different values of Ji/ J = 0.5 and 1.0. As one can see, the critical 
lines of both the lattices terminate at common disorder point as predicted by the ground-state 
analysis. The most obvious difference in the criticality of both IHMs thus rests in a reentrant 
transition of the triangular lattice (Fig.^) observable in the vicinity of disorder points. 

3.2 Model B 

Ground-state behaviour of IHMs with an asymmetric decoration of the triangular (Fig.QJ;) and 
orthogonal-dimer (Fig.QJl) lattices is strongly reminiscent to the one discussed earlier. As a mat- 
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ter of fact, the ground state constitute the aforementioned FP and DP that are bounded by the 
first-order transition fine Ac — \/ {1 + 3Ji/J)(l + Ji/ J). The occurrence of DP can be once 
again related to the spin frustration caused by the Heisenberg spins, which are antiferromag- 
netically aligned in the longitudinal direction and ferromagnetically coupled in the transverse 
direction. In the consequence of that, all Ising spins must be strongly frustrated because they 
interact with the longitudinal components of the Heisenberg spins only. Hence, one may con- 
clude that an appearance of the disordered state |^dp) = Ili ID* Ilfc "T^^II^J-^ IJ-t))* above 
certain Ac (depending on the lattice topology) represents a generic feature of all IHMs having 
the magnetic structure of some close-packed lattice. 

Finally, we shall close our discussion with the finite-temperature phase diagram of IHMs with 
the asymmetric decoration of the triangular and orthogonal-dimer lattices. It is quite interesting 
to ascertain from Fig.|4]that the critical temperatures of both these lattices overlap as they exhibit 
almost the same decline (not discernible within the displayed scale) upon strengthening of A. 

4 Conclusion 

In the present article we have provided an exact solution to a special sub-class of IHMs which 
consist of the Heisenberg spin pairs coupled either to two or three outer Ising spins. The investi- 
gated spin systems exhibit an interesting quantum behaviour manifested in a remarkable kind of 
geometric spin frustration, which appears notwithstanding the purely ferromagnetic interactions 
of the considered model systems. We found a convincing evidence that this behaviour emerges 
due to a competition between the familiar Ising interactions and the easy-plane XXZ Heisenberg 
interactions, which favour ferromagnetic (antiferromagnetic) ordering of the transverse (longitu- 
dinal) components of the Heisenberg spin pairs. The most important question that arises from our 
study is whether this behaviour represents a general feature of all quantum models with mixed 
easy-axis and easy-plane bonds and whether it would persist also in the pure Heisenberg spin 
systems. This question represent a great challenge for our next investigations. 

Acknowledgement: The authors acknowledge financial support provided by Ministry of Educa- 
tion of Slovak Republic given under the grants VEGA 1/2009/05 and APVT 20-005204. 

References 

[1] E. Manousakis: Rev. Mod. Phys. 63 (1991) 1 

[2] S. Sachdev: Quantum Phase Transitions, Cambridge University Press, Cambridge, 1999 
[3] A. P. Ramirez: in Handbook of Magnetic Materials, Vol. 13, 

ed. K. H. J. Buschow, North Holland, Amsterdam, 2001 
[4] J. Strecka, M. Jascur: Phys. Rev. B 66 (2002) 174415 
[5] J. Strecka, M. Jascur: Phys. Stat. Solidi B 233 (2002) R12 
[6] J. Strecka, M. Jascur: J. Magn. Magn. Mater. 272-276 (2004) 987 
[7] I. Syozi: Progr. Theor. Phys. 6 (1951) 341 
[8] L. Onsager: Phys. Rev 65 (1944) 117 
[9] M. E. Fisher: Phys. Rev 113 (1959) 969 
[10] 1. Syozi: in Phase Transitions and Critical Phenomena, Vol. 1, 

eds. C. Domb, M. S. Green; Academic Press, New York, 1972, p.269-329 



